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Abstract 

Using the auxiliary field representation of the simplicial chiral models on a (d - Id- 
dimensional simplex, the simplicial chiral models are generalized through replacing 
the term Tr(AA T ) in the Lagrangian of these models by an arbitrary class function 
of AAt; V{A/tf). This is the same method used in defining the generalized two- 
dimensional Yang-Mills theories (gYM^) from ordinary YM2. We call these models, 
the "generalized simplicial chiral models". Using the results of the one-link integral 
over a U(iV) matrix, the large-iV saddle-point equations for eigenvalue density function 
p(z) in the weak (/? > f3 c ) and strong ((3 < f3 c ) regions are computed. In d = 2, where 
the model is in some sense related to the gYM2 theory, the saddle-point equations are 
solved for p(z) in the two regions, and the explicit value of critical point (3 C is calculated 
for V{B)=TiB n (B = AA*). For V(B)=TtB 2 ,TrB 3 , and TYB 4 , the critical behaviour 
of the model at d = 2 is studied, and by calculating the internal energy, it is shown 
that these models have a third order phase transition. 
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1 Introduction 



One of the useful approaches in studying the large-scale structure of non-Abelian gauge 
theories is via the lattice formulation introduced by Wilson p|, in which a matrix-valued 
field theory is defined on a lattice. After the work of Wilson, several other matrix models 
have been introduced and their properties have been studied. One of the most important 
features of matrix models, is their large-iV behaviour, where N is a quantity related to the 
number of field components. This goes back to the original work by Stanley on the large- 
TV limit of spin system with O(N) symmetry, soon followed by Wilson's suggestion that 
the 1/N expansion may be a valuable alternative in the context of renormalization group 
evaluation of critical exponents, and by 't Hooft extension [[| to gauge theories and, more 
generally, to fields belonging to the adjoint representation of SU(iV) groups. Therefore the 
1/N expansion is probably the most important nonperturbative and analytical tool presently 
available in the study of matrix models. Unfortunately, application of this tool is limited to 
a small number of few-matrix systems. This number is even smaller if we restrict ourselves 
to the case of unitary matrix fields, which is especially relevant to the problem of lattice 
QCD. To the best of our knowledge, the only solved examples are Gross-Witten's single- 
link problem @] and its generalizations, the external field problem ||, and L = 3,4 chiral 
chains [0, §]. 

One of the important classes of the unitary matrix models, as mentioned above, are 
few unitary-matrix systems. The interest for few-matrix models may arise for various rea- 
sons. For example, their large-iV solutions may represent non-trivial bench-marks for new 
methods meant to investigate the large-A" limit of more complex matrix models, such as 
QCD. Furthermore, every matrix system may have a role in the context of two-dimensional 
quantum gravity; indeed, via the double scaling limit, its critical behaviour is connected 
to two-dimensional models of matter coupled to gravity. For more discussion about these 
models see ref. |§. 

An interesting class of finite-lattice matrix models is obtained by considering the pos- 
sibility that each of a finite number of unitary matrices may interact in a fully symmetric 
way with each other, while preserving global chiral invariance; the resulting system can be 
described as a chiral model on a (d — l)-dimensional simplex. A simplex is formed by con- 
necting d vertices by d(d— l)/2 links. These models are known as "simplicial chiral models" 
(SCM) [|l(J. In ref. [11], the large- N saddle-point equations for density function p(z) of 
these models have been found. The main strategy for the determination of these equations, 
is based on the introduction of a single auxiliary variable A (a complex matrix), leading to 
the decoupling of the unitary matrix interaction. The resulting action, contains a Tr^A^) 
term and some other linear terms in A and A^ . Performing the single-link external field in- 



tegral, the authors of |I| have found the saddle-point equations. In d = 2 (where the model 
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corresponds to the Gross- Witten single-link problem, which in turn is equivalent to large- N 
QCD 2 with Wilson action on the lattice), d = 4, and d — > oo, the saddle-point equations 
have been studied analytically. It was shown that the critical value of (3{= {g 2 N)~ l , where 
g is coupling constant) is f3 c — 1/d for all d. Also for < d < 4, it has been shown that 
the model exhibits the third-order phase transition, and for 4 < d < oo it has a first-order 
phase transition. 

On the other hand, it is known that the pure 2-dimensional Yang-Mills theory (YM2), 
with Tr(F 2 ) Lagrangian, can be also represented by the Lagrangian iTi(BF) + Ty(B 2 ) , in 
which B is an auxiliary pseudo-scalar field in the adjoint representation of the gauge group. 
Path integration over B, leaves an effective Lagrangian Tr(F 2 ). Now the generalized 2- 
dimensional Yang-Mills theory (gYM 2 ) is a theory with the above Lagrangian in which the 
term Tr(B 2 ) is replaced by an arbitrary class function f(B) ||12|| . The partition function of 



gYM 2 has been calculated in different contexts in [|L2jl , JT3[, and fl4j| . Also the phase structure 
of these models, at large-iV limit, have been studied in [15[], in |16[ for f(B)=Tr(B 4 ), and in 
T7| for }\B)=Ti(B Q ) and Tr(E 2 ) + #Tr(5 4 ). In pj and 0, it has been shown that these 



models have a third-order phase transition, as ordinary YM 2 . 

In this paper, we want to introduce and study the generalized simplicial chiral models 
(gSCM), in the same fashion that gYM 2 has been defined, that is, through replacement of 
the term Tr(A4''") in SCM by V(AAt'), where V(AAt') is an arbitrary class function of AA'. 
In this way, we have a new few-matrix model on a (d— l)-simplex which, in the special case 
V =Tr(v4v4t), reduces to ordinary SCM. In d — 2, where this model is in some sense related 
to gYM 2 , we see that the study of the phase structure of the model is much easier than the 
corresponding studies in gYM 2 . 

The paper is organized as follows. In section 2 we introduce the gSCM and, using the 
single-link integral method, find the large-iV saddle-point equations for eigenvalue density 
function p(z) in the weak and strong regions. Since these equations for arbitrary V(B) 
(B = AA^) are complicated, they can not be solved analytically for arbitrary d. Note that 
this is also true for V(B) =TrB, i.e. SCM, and as mentioned earlier, for very few cases 
one can compute the density function p(z) analytically. Here the situation is worse and 
therefore we restrict ourselves to d — 2. In section 3, the general case V(B) =Tr(B n ), with 
n an arbitrary positive integer number, is considered and an expression for internal energy 
in terms of p(z) (for arbitrary d), and also the critical value f3 c for d = 2 are found. In 
section 4, considering the leading terms of p(z) (near the critical point) in weak and strong 
regions (again in d — 2), it is shown that for V(B) =Ty(B 2 ), Tr(B 3 ), and Tr(B A ) the model 
exhibits a third order phase transition. We believe that this behaviour (the third order phase 
transition) is the same for all V(B) =Tr(B n ) models. 
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2 gSCM and their saddle— point equation 



If we assign a U(iV) matrix to each vertex of a (d— 1) -dimensional simplex, then the partition 
function of simplicial chiral models is defined by 



Z d (/3,N) = J {[dUiexpiNpfl E Tr(UiU} + U}Uj)}, (1) 

i=l i=l j=i+l 

where dUi is the normalized invariant Haar measure. The d— matrix simplicial model has 
an underlying permutation symmetry instead of the cyclic symmetry of the <i-matrix chiral 
chain. But for d — 1,2, and 3 these two symmetries and the associated models are equivalent. 
The free energy, internal energy, and specific heat are respectively given by 

F d (P,N) = ^lnZ d (P,N), 
11(3 N) _ ^F d (f3,N) 

cMN)=f d -^. (2) 

As mentioned earlier, the main strategy for determination of the large- N saddle-point equa- 
tion is based on the introduction of a single auxiliary variable A (a complex matrix), leading 
to the decoupling of the unitary matrix interaction |10| 

Zd = Zd/ Zq, (3) 

where 

/d 
H dUidAexp{-N(3TiAA^ + N/3TtA ^ U} + N(3TtA^ J2 U i~ N 2 fid}. (4) 

i=l i i 

Performing the single-link integral over the matrices Ui 

e NW(BBi) = J dUexp [ NTr (B^U + U*B)], (5) 

we obtain 

Z d = J dAexp{-N/3TrAA^ + NdW(f3 2 AA^) - N 2 f3d}. (6) 

Now we define the generalized simplicial chiral model (gSCM) through the partition 
function 

Zd,V = Zdy/Z y, (7) 

where 

Z dy = J dAexp{-N(3V{AA^ + NdW{(3 2 AA^) - N 2 (3d}, (8) 
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in which ^(AA^) is an arbitrary class function of AA\ that is 



ViGAA^G- 1 ) = V(AA^), VG e U(N). (9) 
The requirement (^) is satisfied for every polynomial 

V(AA ] ) = a n Tr(AA^) n . (10) 



71=1 



The crucial point in our analysis is that the integrand in eq.@ is a function of eigenvalues 
Xi of the Hermitian semi-positive-definite matrix j3 2 AA^ . Moreover, as Z^v is invariant under 
A — ► G^ AG' (where G and G' are arbitrary unitary matrices), one can perform exactly the 
"angular" integration and reduce the problem to that of integration over iV variables Hl8f . 
Therefore, performing the angular integration, eq.(G|) leads to 

Z d>v = J d^{xi)e^{-N(3V{xi/(3 2 ) + NdWfa) - N 2 f3d}, (11) 

where 

dfi(xi) = Y[dxi Y[(xi - xj) 2 . (12) 

i i>j 

In the large- limit, in which we are interested, the free-energy function W(xi), resulting 
from a single-link integral over a U(iV) matrix, can be extracted by solving the Schwinger- 
Dyson equations. It is written in a simple closed form [||, |(J 



W(xi) = 2 J2 V^7+~c --^ ln(v^T^ + yJ Xj +c)-Nc- -N. (13) 

i ij 

We must distinguish two different phases, a weak-coupling regime where c = and 

— Y^-< 1, (14) 

and a strong-coupling regime where c is dynamically determined by the condition 

4?ot = l (15) 

Therefore, in the large-iV limit the partition function (|TT|) reduces to 

r N 

Z d y= / J[dx z e- S ^^\ (16) 

J 8=1 

with the action 

= -2E ln l^ - x i\ + NPV( Xi /P 2 ) + N 2 (3d 

i>j 
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- Nd{2 E V^+~c - E HV^+~c + ^ + c) - JVc - ^JV}, (17) 

i ij 



and the saddle-point equation dS/dxi = becomes (after multiplying it by {1/N)y/Xi + c 



with the condition x« > 0. To study eq.(|T%D, we use the standard technique which is based 
on the eigenvalue density function. It is however convenient first to introduce a new variable 
Zi with definition 

Zi = ^/x~+~c, (19) 

subject to the condition < \fc < z%. If the eigenvalue variable X{ varies in the interval 
[x a ,Xft], < x a < Xb, then the new variable Zj lies in interval [a, b] where a = y/c + x a , b = 
y/c + Xb, and 

< sfc < a < b. (20) 



In the weak-coupling regime, c = 0, one expects in general a = ^Jx^ > 0, and in the 
strong-coupling regime, one expects x a = so that a = ^fc / Op. 



Now, using the large- N eigenvalue density function p(z), which vanishes outside the in- 
terval [a, b], the saddle-point equation (|T8|) can be turned into the following integral equation 

^ - o- - - = (rb " Irl) ■ (21) 

in which we use the expression ( [370] ) for V and V indicates the principal value of integral. In 
this equation, the parameters a and b must be determined dynamically. The normalization 
condition of p(z) is 

rb 

p(z')dz' = 1. (22) 



In the weak coupling regime, c w = and the condition (13) in the large- N limit becomes 

b dz 'P^l< 2 . (23) 
In the strong coupling regime, where c is 

c s = a 2 , (24) 

the condition ([15]) becomes 

/V^ = 2. (25) 
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Eq. ( pip has a somewhat unconventional form when compared to other integral equations, 
and one must perform a few manipulations to obtain a more familiar equation. To do so, 
we define the function H(z) in the complex-z plane 

H(z) = t ^-dX. (26) 

J a Z — A 

This function is analytic on the entire complex plane except for a cut on the positive real 
axis in the interval [a, b}. Then one has 

H(z±ie) = R(z)Tinp{z) , b < z < a, (27) 

where R(z), from eq . (f2l|) , is 

= "f^^-or^-d-^-Hi-z). (28) 

The constraints (p5|) and (p5|) for weak and strong-coupling regimes, respectively, result in 
H w (0) > -2 and H s (0) = -2. But, from eqs. (g7|) and ©, it can be shown that H(0) must 
satisfy H(0) = -d - ^H(0) - inp{0), or H(0) = -2 - (2iri / d) p(0) . So the weak density 
function must satisfy 

-zp w (0)>0, (29) 
and strong density function must satisfy 

- ^p s (0) = 0. (30) 

Therefore, the analytical properties of H(z) in the two regimes must be such that the con- 
straints (|2D| ) and ([Jl]) are satisfied. Using the standard method of solving the integral 
equations one can show that the expression 

H w (z) = ±y/(z-a)(z-b)i (31) 

has the correct analytical behaviour in weak-coupling region. The contour c in ([JT]) is a 
contour encircling the cut [a, b] and excluding z. Deforming c to a contour around the point 
z and the contour Cqo (a contour at the infinity), one finds 

H w (z) = Ruiz) + ±-.y/(z -a)(z-b)i R ; {X)dX (32) 

2« vv 'JL ( z - X)yJ(X-a)(X-b) 
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Remembering c w = and using eqs.([27|) and (|28|), it can be easily shown that the saddle- 
point equation ( PH) reduces to the following expression for the density function p w (z) 



n,p,q P 

d-2 f b dy p w (y) 



L y + z J(b + y)(y + a) 



} for p>/3 c , (33) 



where 

2m -1) ! , . 

C m = ~ T~ — • 34 

2 m m v ' 



In the strong-coupling regime, one can again show that 
HJz) = 



z / z 



b rdX jX-aR s (X) 



2m V z — a Jc A V A — b z — X 

iX fx_ 

JyJ-bz-x 




(35) 



and, as c s = a 2 (eq.(^4|)), arrive at 




P.(z) = - 1 \l— { E -^l(-^) P r lB m C 2n „ 2m „ 2 zia m b^- m -* 



n,m,p,q 



P2n-1 




-2 



+ ^ f J^.^±±PM } for P < (36) 




in which 

(37) 

with Bq = — 1. It can be seen that in the special case a n = 5 n> i (which corresponds to 
V(A4"I") =Tr(AA^)), the eqs. (|33"D and (j36|) reduce to the corresponding equations in . 



As it is obvious from eqs. (|33"D and fl36|) , here the situation is much more involved than in 



SCM, and the most of analytical calculations done in [0]] can not be done here. One of these 
kinds of calculations is discussed in the next section. Therefore, let us restrict ourselves to 
the case d = 2, for which from eqs. (|33|) and ( |3~6"D the density functions are known. 

Other quantities that must be determined in d = 2 are the values of parameters a and b 



in the both regimes. To find these parameters, we note that at \z\ — > oo, eqs.(|26|) and (|22| 



imply that H(z) — > 1/z or (l/<J(z — a) (z — b) ) H(z) — > 1/z 2 . Therefore, we can expand 

(l/y(z — a)(z — b) ^ (R w (z) — iir p w (z)) (at d — 2), and take the coefficients of l/z and 1/z 2 
equal to and 1, respectively. In this way we find the following equations which must be 
solved to determine a and b in j3 > /3 C 

n 2 n-l C m C2n-m-\O m b 2n ^ m 1 — 2 = 0, (38) 
n,m P 
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and 

ft n 

E -^C m C 2n - m a m b 2n - m - (a + b) = 1. (39) 

n,m P 

In the strong region, as (1/ z)yj (z — a) /(z — b)H(z) — > 1/z 2 at \z\ — > oo, we can again expand 
(1/ ' z)\J(z — a)/(z — b)(R s (z) — mp s {z)) (at d = 2) and take the coefficients of 1/z and 1/z 2 
equal to and 1, respectively, and then find the following equations which determine a and 

b in f3 < (3 C 

E i^l(- a2 ) P ( n ~~ 1 )fi m C 2n _ 2p _ m . 1 a m & 2 "- 2 "- m - 1 + 2 = 0, (40) 

n,m,p P \ P / 

and 

- E iS^(- a2 ) P ( n ~ l )B m C 2n ^ m a m b 2n ~ 2 ^ m + a-b = l. (41) 

n,m,p P \ P J 

Both sets of equations (|38|), (|39| ) and (|40|) , (p|) are too complicated to be solved exactly. 
We discuss about the solutions of these equations in section 4. 



3 Some general results for V =Tr(AA^) 



n 



Let us first derive an expression for the internal energy in terms of p(z) for arbitrary d. If 
we denote the internal energy per unit link by U d y, then the first two equations of (@) lead 
to 

d{d-l) TT ldF dy 1 d ~ ~ 

—2— *V = 2^ = w^T^ v - lnZ °^ (42) 

where d(d — l)/2 is the number of links. Note that ln.Z d y « — S d y{xi}, where Xi S £1X6 the 
solutions of the saddle-point equation (|i~8|) with density functions (|33|) and (|36D in the weak 
and strong regions, respectively, and S d y{xi] is given by eq. (|i"7f) . Therefore, using eq . ([17D 
for V = J2i{xi/(3 2 ) n , the eq.(|2]) reduces to 

2r? — 1 1 8 

d ( d ~ = E^-rf- W2 ^ 0>n . (43) 

Now it is interesting that one can calculate the last term of the above equation exactly. In 
d = 0, eq.fln]) is 



J 0,n 



J® i=l i>j 



i>3 

N 



£ D au ..., aN fdxf' dx N x? ■ ■ • ^e-WELWfl" 



a l''"' a A r ~° 
^\a; = JV(iV-l) 
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N 



y , D ai ... aN I ai ■ ■ ■ I aN , (^4) 

^\ Qi =iV(]V-l) 

where D aii ... j0lN are some unimportant constants and I a is 

dxx a e~ N ^ x/l3 )n 



o 



1 /^n-lX(«+D/n +1 



n \ N n 



T(—). (45) 



Therefore 



Z 0n = ^ 2 (2«-i)/« x (/^-independent terms), (46) 
and from this, eq. fl43|) becomes 

211 — 1 1 1 

d (d - l)U d , n = _ £ a* - d + (- - 2)-. (47) 

In the large- N limit, this equation for weak and strong regimes becomes 

d(d - 1)U$ = I z^p w {z)dz - d+ (- - 2)-, (48) 



and 

d(d - 1)17$ = ^ (^ 2 - a 2 )"^)^ - d + _ 2)^, (49) 

respectively. Again, in the case n — 1, where the model is SCM, the above equations reduce 
to the corresponding equations found in [II|. As can be seen from these relations, the power 
of variables in the integrands makes the integrations difficult to perform, even when the 
density functions are known. Let us show this by an example. It can be shown that in 
d > 4, the critical value of a is different from zero (a c ^ 0) JX1|J , and of course at this point 
Ps = Pw = Pc- Therefore, at the critical point we have 



did - l)[UZ - U%>\ = I [(z 2 - ai) n ~ z 2n ] Pc (z)dz. (50) 



For SCM , where n — 1, the above difference reduces to (using eq.(|22])) 

did- l)p7« - C/i?] = - J = - J ^ 0, (51) 

and this simply proves that there exists a first-order phase transition in SCM with d > 4, 
although we do not know the explicit form of p c (z). But if n is different from 1, the same 
calculation is not possible. For example if n = 2, eq. (|50|) gives 



d{d - 1)[U$ - U<g] = ^K~2 f C Pc (z)z 2 dz\. (52) 
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In this case one must know p c (z) to say something about the result. 

So let us again restrict ourselves to d — 2, and try to find the critical value [3 C for 
V =Tr(AA^) n . As mentioned above, at the critical point p s (z) = p w (z) = pftz) and therefore 
the constraint (^) must also hold for p w (z) at this point, i.e., p™ tlcal (0) = 0. Using (|33|) for 
d = 2 and putting z = 0, we obtain 



Jb e C p C 2n -2 P -2a p c b 2 c n - 2p ~ 2 = 0. 



(53) 



p=0 



As the coefficients of all terms in the above equation are positive, and < a c < b c , the only 
nontrivial solution of (|53|) is 

a c = , for d = 2. (54) 



Therefore, at the critical point eqs.(^0[) and ( |4T| ) reduce to 



n 



jBoC^b 2 ™- 1 + 2 = 



ft 



2n- 



and 



n 



ft 



2/i 



(55) 



(56) 



respectively. From these two equations, one can find the critical values b c and f3 c for 
V =Tr(AA^) n in d = 2 as following 

2n 



c 2n-V 



ftc 



2n 



2n - 1 



n(An - 3)!! 



(57) 



2 2n {2n - I) 

In the case n = 1 (SCM), eq. fl57D gives (3 C = 1/2, which is the correct value of f3 c in the 
Gross- Witten single-link problem. 



4 Phase structure in d=2 

In this section we study the phase structure of some special cases of gSCM in d = 2 with 
V =Tr(AA^) n . 

4.1 ra=l (YM 2 ) 

In the weak-coupling regime, eqs.(^) and (^) are linear equations with solutions a = 
2(3 - ^2ft and 6 = 2/5 + y/2ft, and therefore eq.pfl) leads to 



, „ J(b-z)(z-a) J 2(3 -{z- 2/3) 2 
P { : =1 \z) = ^ £ = V in — ■ (58) 
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In the strong-coupling regime, eqs.([4(]) and (f4T|) result a = 1/2 — (3 and b — 1/2 + /3, and 
therefore eq. (|3~6"D leads to 



P * Uj "7r/3V^-a"7r/5V2^-l + 2/?- 

These equations are the density functions that have been obtained in [[l ij (note that in this 
paper, x^s are the eigenvalues of f3 2 AA\ but in ref. fTT| the eigenvalues of A(3 2 AA^). 



If one computes the internal energies fl4~8|) and (f49|) by using these density functions, it is 
seen that 

Ui n=1) = 1 - ^ , for > i (60) 

and 

UjT 1] =P , for /3 < i (61) 

which are the results obtained in ||. It is obvious from eqs. (|60"D and (|6l|) that in the case 
n — 1, there exists a third-order phase transition. 

4.2 n = 2 gSCM 

In this case the eq.fl5"?D leads to 

c 3' 

fl^>=i(§r- ««> 



First consider the weak-coupling regime. The eqs. (|38|) and ( p9|) give, respectively 

5(a 3 + 6 3 ) + 3ab(a + b) = 16/3 3 , (63) 

and 

35(a 4 + 6 4 ) + 20a6(a 2 + b 2 ) + 18aV = 64/3 3 (a + 6+1). (64) 

These equations can not be solved analytically to obtain a and b, but as we want to study 
the phase structure of the model, it is sufficient to look at the solutions near the critical 
point. Therefore we expand the above equations around a c = and b c = 4/3 up to second 
order, and then find a and b in terms of a = (3 — (3 C , again up to second order. The final 
result is 

a w = -5 2 / 3 2 1 / 3 a--5 1 / 3 2 2 /V + ---, 
4 64 

b w = - + -5 2 / 3 2 1 / 3 a - A 5 i/3 2 2/3 a 2 + . (65) 
3 20 320 V ; 
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Now consider the density function p(z) in the weak-coupling regime for a n = 5 n 2 and d = 2 
(eq.(|D): 



Tip 3 



-(a 2 + b 2 ) + -ab + -(a + b)z + z 2 



yj(b-z)(z-a), 



(66) 



using this function, we calculate (using (£8)) the internal energy (per unit link) in the weak- 
coupling regime 

3 



jj(n=2) 



2/3 4 Ja 

9 



2 12 (3 7 



b 3 

p w {z)z i dz 



4/T 



(67) 



If we substitute the expansions (|65| ) in the above relation, we find the internal energy in 
terms of /3, up to order a 2 , as following 

( n=2) = 279 2/3 1/3 _ 333 5 i/3 2 i/3 + 189 2 (6g) 
400 800 160 V ; 

We can follow the same steps in the strong-coupling regime, this time using eqs.(^), (|41), 
(|36|), and fl4"9"|). The final results are 



as = _I 5 2/3 2 l/3 a _ A 5 l/3 2 2/3 a 2 + . . . 



and therefore 



b s = - + -5 2 / 3 2 x / 3 a + A5 1 / 3 2 1 /3«2 + . 
3 20 80 

5 2/3 2 l/3 + x _ ^ 5 l/3 2 l/3 a + §1 2 + 



279 
400 



300 



40 



c/ („=2)_ f/ ( r =2) = | (/3 _ /3c) 2 + 



(69) 



(70) 



The above equation shows that we have a third-order phase transition for V =Ti(AA^) 2 of 
gSCM in d = 2. 



4.3 n = 3 and n = 4 gSCM 

In the case n = 3, b c = 6/5 and /3 C = (3/80)(16) 4 / 5 (189) 1 / 5 ; and in the case n = 4, b c = 8/7 
and P c = (l/14)3 1 / 7 8 6 / 7 (143) 1 / 7 . The procedure is the same as n = 2 and after some 
calculations, we obtain the final results as following 

^(n=3) _ ^=3) = i^(378) 2 / 5 (/3 - p c f + ■ ■ ■ , (71) 
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16896 V ; KH h ' V ; 

These relations show that for V =Tr(AA^) 3 and Tr(i^) 4 of gSCM in d = 2, the models 
exhibit third-order phase transition, as in the YM theory. We expect that the same behaviour 
exists in all gSCM with V =Tr(AA^) n . 



5 Conclusion 

In this paper the SCM was generalized by a method very similar to the method which have 
been used in the generalization of YM theories. This generalization may be interesting from 
several points of view. First, as mentioned in the introduction, many few matrix-models 
can be solved analytically in the large-A limit. gSCM's are models that can be studied 
analytically in this limit. Second, the gSCM at d = 2 can be treated as the lattice version of 
gYM theory, which is an important two-dimensional candidate of QCD. And third, which 
is somehow related to the second, the study of these models (gSCM) is much more simple 
than the gYM. For instance, the proof of the existence of third-order phase transition in 
f(B) =TyB a gYM theory [[16], which corresponds to our V =Tr(AA' f ) 2 model, is more 



difficult than the proof which is given in section 4.2. 
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